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We use the new method of infinitesimal unitary transformations to calculate zero temperature 
correlation functions in the strong-coupling phase of the anisotropic Kondo model. We find the 
dynamics on all energy scales including the crossover behaviour from weak to strong coupling. The 
integrable structure of the Hamiltonian is not used in our approach. Our method should also be 
useful in other strong-coupling models since few other analytical methods allow the evaluation of 
their correlation functions on all energy scales. 
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I. INTRODUCTION 

The Kondo model, originally introduced by Zener |Q, 
is the most basic Hamiltonian describing the physics of 
dilute magnetic impurities in metals. Despite its appar- 
ent simplicity, it has become a paradigm for complex 
many-body effects, due to its non-trivial strong-coupling 
behaviour at low temperatures. Unprecedented theoret- 
ical activity during the last four decades was triggered 
by the seminal work of Kondo |2| , who showed that per- 
turbation theory for an antiferromagnetic exchange cou- 
pling diverges at low temperatures. Within perturbative 
scaling |^-|| the coupling constant grows continuously 
and eventually diverges. This indicates that the impu- 
rity spin is screened in the ground state. However, the 
diverging coupling constant also implies the breakdown of 
the perturbative scaling approach in the strong-coupling 
phase. Therefore perturbative scaling does not provide a 
systematic and controlled expansion describing the weak 
to strong-coupling crossover. This is a typical situa- 
tion in many strong-coupling problems. The develop- 
ment of new, nonperturbative methods like Wilson's nu- 
merical renormalization group (NRG) U and the Bethe 
ansatz [Q eventually led to a quantitative understand- 
ing of the strong-coupling regime in the Kondo model. 
Despite these successes, however, a simple analytical de- 
scription of the crossover from a free spin to the Fermi 
liquid regime in an RG type framework was still missing. 

In the present work, we will attempt to fill this 
gap by applying Wegner's flow equation method 
which was recently employed successfully to diagonalize 
a related strong-coupling problem, the quantum sine- 
Gordon model H. We will mainly focus on the calcu- 
lation of equilibrium dynamical correlation functions at 
zero temperature as this provides one of the most inter- 
esting new perspectives for this new approach. The cal- 
culation of correlation functions throughout the crossover 
region is notoriously difficult with exact methods build- 
ing on the integrability of the model. There is consider- 
able interest in theoretical tools that allow their deter- 



mination in strong-coupling models, which provides the 
main motivation for our work fl(i|| . Also the integrabil- 
ity of the Kondo model is not used in our method and 
nonintegrable perturbations can be studied as well. 
We consider the anisotropic Kondo Hamiltonian 
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Here S is a spin- 1/2 degree of freedom and 
Sfc^Vv^ are creation and annihilation operators for 
the localized electron state at the spin site. Notice that 
the anisotropic Kondo Hamiltonian is equivalent to the 
dissipative two-state system (spin boson model) |0,12|. 
For the noninteracting conduction band we assume a lin- 
ear dispersion. We will treat the model in its bosonized 
form 13 1 : One introduces the bosonic spin density modes 

a(p) = ^= J2 q ' -fC„ T — el i _i c„i | with the com- 
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mutator [cr(— q), cr(q')] = S qq >L/2Tr for q,q' > 0. Here 
L is the system size and all other nonvanishing com- 
mutators can be derived using a'(q) = <j(— q). Notice 
a(—q)\n >= on the ground state for q > 0. The charge 
density modes in |l]) decouple completely and one only 
has to look at the spin density part of (|l|) |l4| 



H = H 

Here H 
set vp = 



Jil 



V227T 



Ja_ 

Ana 



f e iV2*{0) g- + h 



1. The bosonic field is defined as 



qcr(q) cr(—q). In the following we 
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The parameter a > generates the UV-regularization 
of our model. We are interested in its universal prop- 
erties at energies <C a . The longitudinal spin 
coupling can be eliminated by a unitary transformation 
U = exp (i/xS'^O)) with [i = Ju/y/Tl-K. We arrive at 



1 



H = H +g (V(\o,0)S- + V(-X o ,0)S+) 
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with the coupling constant g = J^/Aira and Ao = \/2 — 
J\\/V22tt. Here vertex operators have been introduced 

V(\,x) = exp (iA$(x)) , 

which are complicated many-body objects for a gen- 
eral value of the scaling dimension A Jl5| . Notice 
that our V(X, x) are not normal-ordered. We will 
also need the Fourier-transformed normalized vertex op- 
erators Ct = a- 1 L- 1 / 2 J dxe ipx V(X,x) with a 2 = 
27ra|pa| A2 -Vr(A 2 ). Notice Cl p |fi) = C p \Q) = for 
p > 0, and (n\C p Cl,\Q) = 8 pp ,0(p), (filCjC^fi) = 
Spp>0(—p) for \ap\ < 1. 

Eq. (|2|) with general initial couplings Ao and go is our 
starting point for the flow equation approach. The uni- 
versal low-energy properties of (Q) are equivalent to those 
of the spin-density part of (Q). As proposed by Wegner 
||, we apply a sequence of infinitesimal unitary trans- 
formations to (Q) in order to make it successively more 
diagonal. One can set this up in a differential formulation 



d B H{B) = [ V (B),H(B)} 
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with an anti-Hermitian generator r/(B). H(B = 0) is our 
initial Hamiltonian (||) and all H(B) as obtained by the 
solution of @ are unitarily equivalent to it. In general 
we will have to truncate our system of equations, so this 
equivalence will only be approximate. In order to gener- 
ate a stable sequence of approximations, we choose rj(B) 
such that first (for small B) offdiagonal matrix elements 
are removed which couple states with large energy differ- 
ences (of order B^ 1 / 2 ), and later more degenerate states. 
This is reminiscent of the energy scale separation under- 
lying perturbative scaling. Similar ideas have indepen- 
dently been developed by Glazek and Wilson [similar- 
ity renormalization scheme) [ fl6| that contain Wegner's 
method as a special case. 

As the flow parameter B grows, the effective spin flip 
interaction in the interaction part H mi (B) will become 
increasingly nonlocal and the scaling dimension of the 
vertex operator will begin to flow. We take this into 
account by writing 

H{B) =H a + J dx g(B, x) (V(X(B), x) S~ + h.c.) (4) 

with g(B,x) — L~ x l 2 ^2 p g p (B) e lpx . Wegner has sug- 
gested the generic choice 17W — [iJ ,-ffint] for the gener- 
ator, as this always removes the off-diagonal matrix ele- 
ments in the above energy-scale separated way ||. Our 
generator has the following more general form 

^ = E^ 1} (cls--c p s+) 
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Here rjp 1 ^ = pa p g p follows from Wegner's canonical 
choice. The second generator part has been introduced 
since it will later allow us to generate transformed Hamil- 
tonians with a particularly simple structure. The precise 

(2) 

form of r] pq will be discussed below. 

In general, new many-particle interactions are gener- 
ated in (^). Our truncation criterion is like in J9| the 
operator product expansion (OPE) for vertex operators. 
For example in the anticommutator of vertex operators 
the OPE leads to Jj|| 

{V(X, x), V(-X, y)} = (1 + * A (x - y) d x $(x) + ...) 

k [1 + i(x - y)/a]» + [1 - i(x - y)la\» ) ' (5) 

We only keep the leading nonvanishing term of this ex- 
pansion in our calculation. This amounts to neglect- 
ing terms with larger scaling dimensions (more irrelevant 
terms) and is therefore reminiscent of an RG approach. 
Higher order terms can be successively taken into account 
in a systematic expansion. However, the present approx- 
imation will already turn out to lead to very accurate 
results. Within our approximation, the primary effect of 
the flow (|) is to generate a new interaction term 

H£l = S z J dxf{x)d x ${x) (6) 

in (|]) with a function f(x) that depends on the param- 
eters g p . (^) is identical to the initial longitudinal spin 
coupling. Our strategy is to eliminate after it is 

generated infinitcsimally by again performing a unitary 
transformation U = exp[ieS z &(0)] with a suitable in- 
finitesimal e. Similar to the initial value of the scaling 
dimension that is determined by the longitudinal Kondo 
coupling, we find a further shift in X(B) due to the elim- 
ination of i/ncwi one derives the following flow equation 

87roA 2 (l-A 2 ) 
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In the sequel we will only look at the strong-coupling 
phase (|J_l| > — Jii for small couplings | Jj_| <c 1) since the 
weak coupling regime is trivial. In the strong-coupling 
phase the flow is always directed towards the Toulouse 
line A = 1 which constitutes the strong-coupling 

fixed point of our approach (see fig. ^). For A = 1 the 
vertex operators in (||) obey fermionic anticommutation 
relations [p~5[ and the Hamiltonian is quadratic (equiva- 
lent to a noninteracting resonant level model [[L7j): The 
OPE series (|^) terminates after the leading term and the 
flow equations close exactly since they can readily di- 
agonalize a quadratic Hamiltonian. We flow to a point 
where our method becomes exact, thereby avoiding the 
usual strong-coupling divergence. This is the key im- 
provement of our approach as compared to perturbative 
scaling, which does not make use of the special properties 
at A = 1. 
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FIG. 1. Flow of the scaling dimension for various cou- 
plings J in the isotropic Kondo model. 

The second new interaction that is generated is a po- 
tential scattering term. Only diagonal terms appear 

b£L(b) = £ M^) (cj c p - c p cl) 



(2*) 

with the following choice of the generator rjpq (B) 
(2) _ \ (p + q) a p a q9p9q 
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mti(B = oo) is a term in order 1/L since oj p (B) = 
0(1/ L) and thus represents the only "impurity effect" in 
the Hamiltonian after elimination of the Kondo coupling 
for B — ► oo in our present approximation 

H(B = oo) = ff + H£l(B = oo) . 

Together with (Q) we thus obtain a closed set of flow 
equations 



d B g P = -P 2 9 P + 2 — V$ 9q 



+ ^9p In ^^-J A9 B A - 2pg p u) p (8) 

d B uJ P =pa 2 p gl . 

Neglecting the momentum dependence induced by ?/ 2 ) 
and the potential scattering terms lo p (which become 
important only at energies of the order of the Kondo 
scale), we can parametrize the solution of (||) as g P (B) = 
g{B)e- p B with the running coupling g(B). g(B) is in 
general strongly renormalized. In the low energy limit the 
scaling dimension approaches A = 1 and we can then de- 
duce the Kondo scale from the effective Toulouse Hamil- 
tonian where Tk ~ g(B = oo) 2 /a. For the symmetric 
case J|| = J± = J and a small antiferromagnetic cou- 
pling < J < 1 one e.g. shows T K ~ ^ exp (-^f) , 
where r = (1 + 7 - ln2)/3 w 0.295. NRG and third or- 
der scaling results by Wilson predict the same leading 
exponential behaviour but a slightly different prefactor 
with r = 1/2 pi. 



Our main focus in this work is the calculation 
of the dynamical impurity susceptibility \(t) = 
i6(t) ([S z (t),S z (0)]) at zero temperature. The key ideas 
for the calculation of correlation functions for impurity 
systems within the flow equation approach have been de- 
veloped in Jl9| |: we need to evaluate matrix elements of 
the observable O — S z with respect to the eigenstates 
of H(B = 00), which is simplified by the fact that the 
final Hamiltonian is diagonal. In order to make use of 
this simplification, we need to perform the same series of 
infinitesimal unitary transformations on O as on H(B) 



d B 0(B) = [ V (B),0(B)} 



(9) 



Usually, keeping the Hamiltonian flow (|3|) simple leads to 
a complicated structure of the transformed observable. 
In the present case, however, it turns out that the OPE 
can be used as a unifying truncation criterion for both 
H(B) and 0(B). Using S z = [S+,S~}/2 and the ansatz 

S+(B) = h(B) S+ + S Z Y d p (B) Cl 



for the flow of the observable [M another closed set of 
equations is obtained: 



d B h = - p a p dp 
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d B d p = ~2pg p a p h + 2 Vpq d q 
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We have numerically solved this system of equations in 
combination with (||). It is known from the theory of 
dissipative quantum systems |l9[ that in the thermo- 
dynamic limit the impurity observable S + decays com- 
pletely, i.e. h(B = 00) = 0. With S Z (B = 00) com- 
pletely expressed in terms of band operators, it is then 

trivial to calculate the dynamical susceptibility: since 
( 2^ 

Hn CW (B = 00) is of order 1/L, it does not contribute to 
the dynamics and we need only study the time evolution 
of S Z (B — 00) under Hq. The imaginary part of the 
Fourier-transform of x(t) follows easily 
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Results are shown in fig. |^: The curves display a broad 
maximum at an energy of the order of the Kondo scale 
and a power law decay x"( w ) ~ tj~( 3 ~ 2Q ) at high frequen- 
cies, consistent with results obtained for the spin boson 
model |plf . At low frequencies one finds x"( w ) ~ w > 
where the corresponding slope scales as T^ 2 with a pref- 
actor of order one [p2[ . The normalization condition 
/ °° duj x"( w ) = 7r/4 is fulfilled due to the sum rule 
h2 ( B ) + J2p d l( B )/^ = 1 that holds exactly @. 
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FIG. 2. Universal scaling forms of the dynamical im- 
purity susceptibility for different dissipation strengths 
a = (1 — J||/47r) 2 in the limit of small coupling J±: u>o is 
defined by the maximum of the curves. One finds loq oc Tk- 

The real part can be obtained by a Kramers- 

Kronig transformation. In particular, the static suscep- 
tibility in response to a local field is given by \o = x'( w = 
0) and therefore contains spectral information on all en- 
ergy scales. In the limit of small spin flip coupling J± 
we have numerically verified the power law behaviour 
Xo ~ ^ known from the ohmic spin boson model 

with the dissipation parameter a = (I — J^/Att) 2 

Summing up, we have applied Wegner's flow equations 
to the anisotropic Kondo problem. The Hamiltonian 
is written in terms of vertex operators with scaling di- 
mensions that flow to the Toulouse point. In contrast 
to the perturbative scaling approach, our flow equations 
close exactly at the Toulouse point with finite couplings 
and a nontrivial strong-coupling behaviour. Similar to 
the analysis of the sine-Gordon model || no strong- 
coupling divergence of the coupling constants is encoun- 
tered. We are thus able to give an analytic description of 
the crossover from weak to strong coupling in a system- 
atic expansion that can be improved by taking higher 
orders in the OPE into account. The stability of the 
strong-coupling fixed point cannot be endangered in this 
expansion. 

In our solution of the flow equations we have focussed 
on the equilibrium dynamics of the impurity at zero tem- 
perature. As an example we have calculated the universal 
scaling forms of the local dynamical susceptibility x(w) 
and found good agreement with results known from the 
spin-boson model. This exemplifies the usefulness of our 
approach for the calculation of correlation functions in 
strong-coupling problems, and should be of interest in 
other models as well. Notice that we have not used the 
integrable structure of the Kondo model in our method. 
Future work will concentrate on the effect of finite tem- 
perature, and on the calculation of static impurity contri- 
butions where the newly generated potential scattering 
terms give the leading contribution (in contrast to the 
dynamical susceptibility where these terms play no role) 



and our present approximation must be improved. 
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